Chandra's Discovery of Activity in the Quiescent Nuclear Black 

Hole of NGC 821 



G. Fabbiano 1 , A. Baldi 1 , S. Pellegrini 2 , A. Siemiginowska 1 , M. Elvis 1 , A. Zezas 1 , J. 

McDowell 1 

1 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138; 
gfabbiano@cfa.harvard.edu; abaldi@cfa.harvard.edu, aneta@cfa.harvard.edu, 
elvis@cf a. harvard, edu, azezas@cfa. harvard, edu, jcm@cfa. harvard, edu 

2 Dipartimento di Astronomia, Universita' di Bologna, via Ranzani 1, 40127 Bologna 

( Italy ); silvia.pellegrini@unibo.it 

May 18, 2004 

ABSTRACT 

We report the results of the Chandra ACIS-S observations of the elliptical 
galaxy NGC 821, which harbors a supermassive nuclear black hole (of 3.5 x 
10 7 M Q ), but does not show sign of AGN activity. A small, 8.5" long (~ 1 kpc 
at the galaxy's distance of 23 Mpc), S-shaped, jet-like feature centered on the 
nucleus is detected in the 38 ksec ACIS-S integrated exposure of this region. 
The luminosity of this feature is Lx ~ 2.6 x 10 39 ergs s _1 (0.3-10 keV), and 
its spectrum is hard (described by a power- law of Y = 1.8tol; or by thermal 
emission with kT > 2 keV). We discuss two possibilities for the origin of this 
feature: (1) a low-luminosity X-ray jet, or (2) a hot shocked gas. In either case, 
it is a clear indication of nuclear activity, detectable only in the X-ray band. 
Steady spherical accretion of the mass losses from the central stellar cusp within 
the accretion radius, when coupled to a high radiative efficiency, already provides 
a power source exceeding the observed radiative losses from the nuclear region. 

Subject headings: GALAXIES: NGC 821 - GALAXIES: NUCLEI - X-RAY: GALAX- 
IES 



1. Introduction 



Luminous quasars, radio galaxies, and Seyfert galaxies have long been associated with 
accretion onto a massive black hole (see review Rees 1984), and the lack of nuclear emission 
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in most galaxies was alternatively debated as evidence for the lack of such a nuclear black 
hole, or for the lack of fuel reaching the hole (Phinney 1983). 

We now know that virtually all galaxies host supermassive black holes (SMBH) in 
their nuclei. High resolution observations of the nuclei of elliptical galaxies and bulges have 
established the presence of these SMBHs (Richstone et al 1998; Magorrian et al 1998; van der 
Marel 1999), whose mass is loosely correlated with the galaxy/bulge luminosity (Magorrian 
et al 1998) and tightly correlated with the central velocity dispersion (Ferrarese & Merritt 
2000; Gebhardt et al 2000) and the central light concentration (Graham et al 2001). These 
results remove the black hole variable from the equation, and make the absence of luminous 
AGN in these galaxies more puzzling. Although a few faint AGNs have been detected in 
X-rays (e.g, in the L X /L E ~ 10~( 6 ~ 7 ) range, where L E is the Eddington luminosity of the 
SMBH; Fabbiano & Juda 1997; Pellegrini et al. 2003a; Ho et al. 2001; Fabbiano et al. 2003), 
we still do not have a clear picture of what impedes the formation of a luminous AGN, and 
of how this is related to the existence and physics of a circum-nuclear accretion disk (e.g., 
Martini et al. 2001). Given that elliptical galaxies tend to host large quantities of centrally 
concentrated hot interstellar medium (ISM; see Fabbiano 1989 for an early review), lack of 
fuel does not seem to be an option. 

For example, the circum-nuclear regions studied so far with the Chandra ACIS show the 
presence of hot gas close to the accretion radius; this implies, when applying the spherical 
Bondi accretion theory, accretion luminosities values comparable to those of luminous AGNs, 
if the gas close to the SMBH joins an accretion disc with a standard radiative efficiency of 
~ 0.1 (e.g., Loewenstein et al. 2001; Fabbiano et al. 2003). Since these luminosities are 
not observed, the options left are that the radiative efficiency is orders of magnitude lower, 
as in an advection dominated accretion flow (ADAF) and its variants (e.g., Narayan 2002), 
or that accretion is unsteady (Binney & Tabor 1995, Ciotti & Ostriker 2001). In addition, 
nuclear jets can carry away a large fraction of the estimated accretion power, a possibility 
that has been found very reasonable for M87 (Di Matteo et al. 2003), IC4296 (Pellegrini et 
al. 2003b) and IC1459 (Fabbiano et al. 2003). NGC 821 represents another promising case 
to test the hypotheses above. 

NGC 821 is an E6 galaxy with possible disky isophotes and centrally peaked (power- 
law) surface brightness (Ravindranath et al 2001). At a distance of 23 Mpc (e.g., de Zeeuw 
et al. 2002), the ACIS resolution at the aim point corresponds to 55 pc. Prior to our 
Chandra observation, NGC 821 had been observed, but not detected, in X-rays with ROSAT 
(< 5 x 10 40 erg s -1 ; Beuing et al. 1999). This limit is ~ 10 5 times below the Eddington 
luminosity of the nucleus, based on the SMBH mass of 2.8—5.8 x 1O 7 M , which was measured 
from HST (STIS) spectra of the nuclear region coupled to dynamical modeling (Gebhardt 
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et al. 2003). NGC 821 was observed with Chandra as part of a mini-sample of extremely 
faint SMBHs extracted from the list of Ho (2002), which will be the subject of a follow-up 
paper. As reported in Ho (2002) and Ho et al. (2003), NGC 821 has not been detected 
in radio continuum nor in optical emission lines (Ha and H/3), and is a good example of 
quiescent SMBH. Here we report the results of the Chandra ACIS observation of NGC 821, 
that has led to the discovery of an S-shaped feature, suggestive of either a weak two-sided 
X-ray nuclear jet, or of hot shocked gas. 

2. Observations and Data Analysis 

NGC 821 was observed with Chandra (Weisskopf et al 2000) ACIS-S (PI: Fabbiano) 
on November 26, 2002 (ObsID: 4408) and on December 1, 2002 (ObsID: 4006) for a total 
exposure time of 38 ks. Table 1 is a summary of the relevant properties of NGC 821 and of 
the observing log. 

No significant background flares were observed in these data, so no further screening was 
necessary. A time-dependent gain correction (URL: http://cxc.harvard.edu/contrib/alexey/tgain/tgain.html) 
was applied to the Standard Data Processing Level 2 event files, before further analysis. The 
data were then analyzed using the CIAO v3.0.1 software (CALDB 2.23) and XSPEC vll.2.0 
for the spectral fits. 

2.1. X-ray image 

The two observations were merged, taking into account the relative aspect solution of 
the two data sets, with the CIAO task mergejxll. From the resulting dataset, images were 
extracted in three spectral bands (Red = 0.5-1 keV; Green =1-2 keV; Blue = 2-4 keV), 
covering the spectral range in which most source counts are detected. The images were then 



Table 1. NGC 821: Properties' 1 and Chandra ACIS-S Observation Log 



m° Bt 

(mag) 
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M. 
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Texp. 

(ks) 


-20.71 


23 


2.6 
209 d 


6.2xl0 2() 


<5xl0 40 b 
< 1.3 x 10 3Se 


2.8-5.8 c 


55 


4408 
4006 


Nov. 26, 2002 
Dec. 1, 2002 


24.6 
13.3 



a Unless otherwise noted, the galaxy properties are as listed in NED (NASA / IPAC Extragalactic Database); b Beuing et al 



1999; c Gebhardt et al. 2003, rescaled for the distance adopted here; d Prugnicl & Simien 1996; e Ho et al. 2003 
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smoothed using csmooth with scales ranging from 1 to 20 pixels (0.5"- 10"). Fig. 1 shows the 
resulting 'true-color' X-ray image of a 3' field centered on NGC 821. Given the small region 
considered, no exposure correction was needed. The ellipse indicates the -D25 isophote from 
RC3. Most of the X-ray emission is concentrated in the central regions of the galaxy, and it 
includes a diffuse component, a few point-like sources, and a brighter central/nuclear region. 

A higher resolution 'true-color' image of the central region of NGC 821 is shown in 
fig. 2a, where the data are displayed at the original observed resolution, without smoothing. 
This figure shows clearly a north-south elongated, hard central feature. Within the Chan- 
dra aspect uncertainties (< 0.5"), this feature is centered on the nucleus of NGC 821, at 
RA=02 /l 08 m 21.14 s , Dec=+10°59'41.7" (J2000, with uncertainty of 1.25"; from the 2MASS 
survey, as reported in NED). While a point-like source may be superposed on the southern tip 
of this elongated feature, the general form is suggestive of a two-sided bent jet, or S-shaped 
filament centered on the nucleus. This feature is approximately 8.5" long, corresponding to 
~ 1 kpc at the distance of NGC 821. 

2.2. Point source detection and spatial analysis of the nuclear feature 

We ran the CIAO wavdetect tool on the full-band ACIS image, with scales between 2 
and 4 pixels (1-2"), and detected 11 sources (above 3 a) in the area shown in fig. 1. We 
assume a power-law spectrum with T = 1.8 and Galactic N H = 6.4 x 10 20 cm~ 2 (Stark et al 
1992), consistent with the emission of low-mass X-ray binaries (LMXBs) detected in early- 
type galaxies (see Kim & Fabbiano 2003). With this spectrum, the detected sources have 
0.3-10 keV 'emitted' fluxes in the 1.4 x 10 -15 — 1.4 x 10 -14 erg cm~ 2 s _1 range, corresponding 
to luminosities of 1.2 x 10 38 — 1.2 x 10 39 erg s _1 , if they indeed belong to NGC 821. The 
faintest sources we can detect in our exposure of NGC 821 are at the upper end of the 
luminosity distribution of the populations of LMXBs detected in other elliptical galaxies 
with Chandra (see e.g., Kim & Fabbiano 2004). Based on these results, we expect that most 
LMXBs in NGC 821 will have L x < 1 x 10 38 erg s" 1 and therefore will be undetectable. 
These LMXBs, however, will contribute to the unresolved 'diffuse' galaxian emission. 

In the central region shown in fig. 2a, there are four sources: the isolated point-like 
source at the north-east of the central complex (source NE, with 0.3-10 keV L x = 4.7 x 10 38 
erg s _1 ), and three sources in the central elongated emission region; these are identified by 
ellipses in fig. 2b, and named SI, S2 and S3, going from north to south. Further analysis 
(see below) demonstrates that these three central emission regions are not point-like and 
therefore cannot be explained with the serendipitous positioning of three luminous galaxian 
LMXBs in NGC 821. No discernible features can be seen in the optical image of NGC 821 at 



- 5 - 



the positions of these three sources. Using the Deep Survey source counts in the 0.5-2 keV 
band (Rosati et al 2002), the number of expected sources at fluxes of > 1.4 x 10~ 15 erg cm -2 
s" 1 is < 8 x 10~ 3 , so that chance detection of background sources is very unlikely. We cannot 
exclude a peculiar clustering of very luminous galaxian X-ray binaries, and future deeper 
data will be needed to explore this point further. 

To establish the spatial properties of sources SI, S2 and S3, we have compared the 
spatial distribution of counts from each of them with that of the on-axis image of the quasar 
GB 1508+5714 (Siemiginowska et al 2003a). GB 1508+5714 can be used as a good repre- 
sentation of the Chandra ACIS-S PSF for our analysis, since the sources we are interested 
in are also at the aim point and similarly hard. With a count rate of ~ 0.05 count s _1 , 
the image of GB 1508+5714 (ObsID 2241) is not affected by pile-up, and contains 5,300 
counts within 2" of the centroid of the count distribution. From this image, we deter- 
mine the ratio of counts within the l"-2" annulus to those in the central 1" radius circle 
to be Ratio(PSF)= 0.043 ± 0.001 (la). The isolated NE source in this central field yields 
Ratio(NE)= 0.057 ± 0.054 (from a total of 36 source counts), entirely consistent with that 
of our reference quasar, confirming that GB 1508+5714 gives a good representation of the 
PSF. Instead, the analogous ratios for the background-subtracted counts from SI, S2 and 
S3 demonstrate that the emission is extended in all cases. Using the wavdetect centroids, 
we obtain Ratio(Sl)= 0.81 ± 0.26, Ratio{S2)= 0.94 ± 0.27, and RaUo{S3)= 0.29 ± 0.10. The 
total number of source counts in the three cases are 56, 67, and 71, respectively, larger than 
for source NE. This comparison demonstrates that the spatial distributions of the source 
counts from SI and S2 are definitely not consistent with the PSF. Although the spatial 
count distribution of S3 is more peaked and could contain a point-like component, diffuse 
emission is also present. This analysis reinforces our suggestion that the central emission 
feature is intrinsically elongated. 

Given that the central emission (S2) is not consistent with a point-source, we can only 
estimate a 3 o upper limit on the luminosity of a nuclear point-like AGN. We used a 1" x 1" 
(2x2 pixels) sliding cell over the entire area covered by the S-shaped feature, and assumed 
as the background level the maximum value detected in the sliding cell. This conservative 
approach yields a 0.3-10 keV L x < 4.2 x 10 38 erg s^ 1 for a power-law spectrum with T = 1.8 
and line of sight Galactic N H . This limit indicates that a central point-like AGN would have 
a luminosity not exceeding that of normal LMXBs. 
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2.3. Spectral analysis 

We extracted spectral counts both from the S-shaped feature (regions in fig. 2b) and 
the surrounding diffuse emission. The latter was taken from a 20" circular region centered 
on the nucleus, excluding the S-shaped feature and the NE source. The field background 
was extracted from a source-free 50" radius circular area 1.8' from the nucleus of NGC 821. 
Spectral counts were extracted separately for the two observations using the CIAO script 
acisspec, which extracts a spectrum for both the source and the background region and cre- 
ates weighted response matrices and ancillary response files (ARF). The ARFs were corrected 
for the time-dependent degradation in the ACIS quantum efficiency, using the CIAO script 
apply _acisabs. The source spectra from the two observations were then coadded (using the 
FTOOL mathphd) to optimize the signal-to-noise ratio of the data. We used the FTOOLS 
addarf and addrmf to combine the responses with their appropriate weights. 

Hardness ratios (HR1=M-S/M+S; HR2=H-M/H+M; where S=0.5-l keV, M=l-2 keV, 
and H=2-4 keV) are plotted (with 1 a statistical errors) in fig. 3, and compared with power- 
law and Raymond-Smith emission models. The hardness ratios of the SI, S2, S3 regions 
are all consistent with hard emission, either a power-law spectrum with T between ~ 1 and 
~ 2.2, or thermal emission with kT ~ 2 — 20 keV. The diffuse X-ray emission of an elliptical 
galaxy is the combination of the soft emission of the hot ISM and of the hard emission of 
the population of LMXBs below our detection threshold (see e.g, Kim & Fabbiano 2003). In 
NGC 821 the hardness ratios of the diffuse emission are hard, suggesting a dominant LMXB 
component and little hot ISM. 

With XSPEC, we fitted the data in the 0.3-10 keV energy range, rebinned to have at 
least 15 counts in each energy bin. For the S-shaped feature we adopted an absorbed power- 
law model (XSPEC model: wabs(wabs(pow))), with N H consisting of both a Galactic and 
an intrinsic component. The results of the fits are listed in Table 2, with 90% errors on one 
significant parameter. The S-shaped emission is well fitted with a power-law spectrum with 
T ~ 1.8, typical of AGN spectra, although the uncertainties are large. As discussed in § 2.2, 
there could be a possibly unrelated point-like source at the southern tip of the S-shaped 
feature. Given the very few counts detected in our observations, we were unable to perform 
a meaningful spectral analysis, eliminating this source. However, the X-ray hardness ratios 
of fig. 3 show that the X-ray colors for the three regions SI, S2, S3 are remarkably similar, 
and consistent with the power-law model derived from the minimum-x 2 fit. 

For the diffuse emission, following e.g. Kim & Fabbiano (2003), an optically thin thermal 
component was added to the power-law model (XSPEC model: wabs(wabs(apec+pow))). The 
metal abundance of the thermal component was fixed to 0.3 times the solar value (keeping 
the elemental solar ratios of Anders & Grevesse 1989). The spectrum of the diffuse galaxian 
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emission (see Table 2) is consistent with a hot gas with a temperature kT ~ 0.5 keV, or 
cooler, typical of X-ray-faint early-type galaxy halos (e.g. Pellegrini & Fabbiano 1994, Irwin 
& Sarazin 1998). The power-law component, although ill-defined, is needed to obtain an 
acceptable fit for the diffuse emission, as clearly suggested by fig. 3; this is in agreement 
with the presence of an unresolved LMXB population. The x 2 f° r a simple thermal model 
(APEC; Smith et al 2001) is 25.8 with 12 degrees of freedom. Comparing this with the 
two-component x 2 by means of an F-test, we obtain a chance probability of 10 -3 . Fitting 
the diffuse emission with a single power-law model, also results in a worse x 2 ■, but in this 
case the F-test probability is only 12%. While the two-component model is plausible for the 
diffuse emission, the rather low signal-to-noise ratio of our data does not allow a stronger 
discrimination among models. In the following discussion we will use the results of the two- 
component model fit; however, future deeper observation of NGC 821 are essential to firmly 
constrain the characteristics of the hot ISM in this galaxy. 

Table 2 also lists the best-fit unabsorbed luminosities for the S-shaped feature, the total 
diffuse emission and the thermal component of the diffuse emission. The latter is only 10% 
of the total diffuse emission, indicating that NGC 821 is singularly devoided of hot ISM. 
Strictly speaking, we have a (90%) upper limit of 3.5xl0 38 erg s" 1 on the luminosity of the 
gaseous component. The uncertainties on these luminosities reflect the uncertainties on the 
spectral parameters. From the emission measure of the soft thermal component we estimate 
an 'indicative' electron density n = 4.1^° 5 9 x 10~ 3 cm~ 3 . While future deeper data are needed 
for a definite measure, the above estimate is useful for gaining first-cut insights on the nature 
of the S-shape feature, and in this spirit we will use it in the following discussion. 



Table 2. Results of Spectral Fits 



Net Counts 2 N H kT Lx(0.3 - 10 keV) 
(0.3-10 keV) X 1 0/ (xlO 21 cm' 2 ) (keV) (erg s' 1 ) 

S-shaped 141 6.7/5 1.41+^ • • • 1.82±g's6 2 - 6 ± °- 5 x 1q39 

6.5/5 0.75^ 5.14^00 ... " 1.9^ x 10 39 

Diffuse 174 6.5/10 < 9.00 0.46^;| 1.27lJ;J| 3.5 ± 0.6 x 10 39 

Thermal ••• ••• ••• ••• ••• 3.4±0.1xl0 38 
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3. Discussion 

NGC 821 was observed with Chandra because the existing ROSAT upper limit on the 
X-ray luminosity of its nuclear SMBH implied extremely sub-Eddington nuclear emission. 
Our observations fail to detect a point-like source at the nucleus, down to a 3a limit of 
(0.3- 10 keV) < 4.2xl0 38 erg s" 1 , ~ 100 times fainter than the ROSAT limit (see Table 1). 
The nucleus is not detected in the FIR or in H 2 (Georgakaki et al. 2001) arguing against 
a strongly obscured AGN (and against a nuclear starburst). There is also no sign of strong 
intrinsic absorbing column in the X-rays (see Table 2). The general AGN 'quiescent state' is 
also supported by the lack of radio continuum and optical line emission (Ho 2002; Ho et al. 
2003). We detect instead an elongated, possibly bent, emission feature, strongly suggestive 
of a two-sided X-ray jet or S-shaped filament. This feature, extending over ~ 1 kpc, has a 
hard spectrum consistent with a power-law with best-fit T ~ l-8^o;<3 (a—O.S) or, if thermal, 
kT > 2 keV. The X-ray luminosity of this feature is ~ 1.9 — 2.6 x 10 39 ergs s" 1 , corresponding 
to ~ 5 x 10~ 7 of the Eddington luminosity of the SMBH. 

The nature of this intriguing elongated feature is investigated below, considering the 
possibilities of a two-sided jet (§ 3.1) and of thermal emission from gas heated by some form 
of energy deposition resulting from nuclear activity (§ 3.2 and 3.3). 

3.1. Is a jet directly emitting the hard X-rays? 

We discuss here the possibility that the S-shaped, hard emission centered on the nu- 
cleus of NGC 821 is a two-sided nuclear jet, as in radio galaxies. The spectral power-law 
slope of this emission has large uncertainties (Table 2), but is consistent with the X-ray 
spectra of other jets (Siemiginowska et al. 2003b; Sambruna et al 2004). However, un- 
like other extragalactic X-ray jets seen in luminous AGNs, where L x (jet)/L X (AGN) ~ 
1-15% (e.g., Schwartz et al. 2003), this 'jet' has no associated core X-ray source, implying 
Lx(jet) / L X (AGN) > 6. This lack of core emission also differs from other X-ray weak AGNs, 
where only a point X-ray source has been seen, typically in nuclei with radio detections, and 
these are usually modelled as Comptonized jets (Baganoff et al. 2001, Markoff et al. 2001, 
Fabbiano et al 2003). An exception is M 87, with its core-jet structure (Di Matteo, Carilli 
& Fabian 2001). The NGC 821 'jet' could be similar to the somewhat steeper spectrum 
(r ~2.3) M87 jet (Wilson & Yang 2002), where the nuclear point-like AGN is fainter com- 
pared with the jet (L x (jet) / L x (AGN) ~ 2). The 0.5 kpc (one-sided) scale of the NGC 821 
'jet' is also similar to the ~1.5 kpc of the M 87 jet, while most radio/X-ray jets extend over 
much larger distances, up to 300 kpc (Siemiginowska et al. 2002, 2003b). The NGC 821 'jet' 
is two-sided, putting a limit on relativistic beaming and suggesting only a relatively slow 
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expansion (v < 0.3c; see e.g. Leahy 1991). The typical expansion velocity of the M87 jet is 
0.5c, though small features can reach 4-6c (Biretta et al. 1995, 1999). 

If a jet is present, we can make some simple estimates of the activity timescale and of 
the energetics involved. First, based on the above limit on the expansion velocity (v < 0.3c), 
we can relate the size of the jet to the activity timescale, finding that the jet may have been 
propagating for a few thousand years, assuming free unimpeded propagation. However, it 
is more likely that the jet is disrupted by interaction with the surrounding hot ISM (n = 
4.ll^° 5 9 x 10~ 3 cm~ 3 ) as, for example, may be occuring in NGC 1316, which has a similar ISM 
density (0.01 cm -3 in the jet disruption region; Kim, Fabbiano & Mackie 1998). Theoretical 
simulations (De Young 1993) of jet propagation into a homogeneous ISM indicate that hot 
ISM with densities of 0.01 particles cm -3 , similar to those of NGC 821, can significantly slow 
down a jet 10 5 times more powerful than the one we may have in NGC 821. 

Considering that the jet is likely to interact with the surrounding hot ISM, its size 
may be indicative of equilibrium between the ram pressure of the expanding jet and the 
thermal pressure of the hot ISM. We estimate the thermal pressure of the ISM from the 
diffuse X-ray emission to be ~ 1.5 x 10~ n dyn cur 2 (P th =2n e kT). Equating this pressure 
to the ram pressure P ram = Pismv 2 , where Pism is the mass density of the ISM, we derive 
a jet expansion velocity v of 0.043c, consistent with the upper limit mentioned earlier, and 
corresponding to a jet expansion (or activity) timescale of 3.8 xlO 4 years. The jet kinetic 
energy can be estimated from E^ in = 0.5P ram V, where V is the volume occupied by the 
jet. To calculate V, we assumed a jet thickness of 10 pc and a length of 1 kpc, obtaining 
V =2.7xl0 60 cm 3 . The kinetic energy input to this volume is then 4.1 xlO 49 erg, and by 
dividing this E^ n by the jet expansion time we derive a kinetic luminosity of 1.7 x 10 37 
erg s _1 . This is ~ 100 times smaller than the X-ray luminosity. However, this may be a 
lower limit if the jet is "overpressured" , as in M87 where the radio cavities in the cluster gas 
suggest overpressurization by a factor ~ 3 with respect to the cluster gas (Reynolds 1997; 
Reynolds et al 1996). However, M87 has a jet kinetic luminosity of ~ 10 43 erg s -1 , and a 
radiative power of ~ 10 42 erg s" 1 , far more powerful jet than in NGC 821. 

What would be the emission mechanism responsible for the S-shaped emission of NGC 821, 
in the jet hypothesis? Synchrotron emission is plausible: a rac iio-x < 0.7 (estimated from the 
radio flux limit of 0.5 mJy at 5 GHz and the X-ray flux at 1 keV); this a ra dio-x is consistent 
with the X-ray slope, and is typical of the synchrotron slope observed in radio lobes (Peter- 
son 1997). The usual problem of short synchrotron lifetime for X-ray producing electrons 
(Feigelson et al. 1981), imply local particle re-acceleration 0.5 kpc from the nucleus. A radio 
detection not far below the current limit is expected if synchrotron radiation produces the 
S-shaped emission at the center of NGC 821. Also in this respect, the NGC 821 emission 
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could be consistent with the M 87 jet: Wilson & Yang find a ra dio-x ~ 0.9 for the knots of 
the M87 jet, and they prefer a synchrotron origin for the X-ray emission of these knots. 

However a ra dio-x < 0.7 is also consistent with the values of a ra( ao-x— 0-7-0.8 reported 
for knots in powerful jets by Sambruna et al. (2004) who favor an Inverse Compton origin for 
the X-ray flux for most of the jet knots, based on the X-ray fluxes lying above an extrapolation 
of the radio-optical slope (a_Ro > «ox)- So the value of a ra dio-x < 0.7 may be merely a 
coincidence, and is in any case only a limit, so it is worthwhile to consider an Inverse Compton 
origin. The seed photons for Comptonization could be either from the synchrotron photons 
within the jet (the 'self-Compton' case) or could be the external to the jet photon field. In 
the self-Compton process the ratio of the synchrotron to Inverse Compton luminosities is 
given by the ratio of energy densities of the magnetic to the synchrotron radiation field. In 
equipartition both luminosities are of the same order, and since the observed X-ray luminosity 
is a factor of at least 10 3 higher than the radio upper limit, the self-Compton case is excluded. 

External Comptonization is also ruled out: Felten & Morrison (1966, eq.47) showed that 
Is/Ic — U B /Up h (i>c I ' v s ) { ^~ rri) l 2 . Here Jg and I c are intensities of synchrotron and inverse 
Compton emission, Ub is the energy density of magnetic field, U p h the energy density of the 
external photon field, v c jv s is the ratio between the frequency of the Compton scattered 
photon and the frequency of the synchrotron photon, m is the power law index of the 
electron distribution, which is linked to the spectral index, a s of the synchrotron emission so 
m — 1 — 2a s . We measure Ic as the X-ray luminosity, and U p h from the starlight from the 
central cusp of NGC 821 (Gebhardt et al. 2003, rescaling their V-band luminosity density for 
the distance in Table 1). The maximum synchrotron luminosity, Is comes from the 0.5 mJy 
flux limit at 5 GHz, integrating over the 10 7 — 10 9 GHz range with a slope of 0.7, and it 
is <4.5xl0 36 erg s _1 . Ub can be derived from the radio upper limit assuming equipartition 
of the magnetic field and the electrons. Because the jet expansion is not relativistic we can 
assume the bulk jet Lorentz factor r ~ 1, simplifying comparison of the radiation fields. 
Assuming equipartition (Burbidge 1959) between the Ub and the electron energy density, 
gives the minimum magnetic field (Krolik, 1999, eq. 9.21), which can be expressed as: 
B min = 4819.84(Lo/V) 2/ V 1/7 (^)) (1 ~ 2aa)/7 G, wher e in our case v = 5 GHz, v x is the lowest 
frequency of the synchrotron emission, Lo is the monochromatic luminosity at 5 GHz, and 
V is the emission volume. We obtain: B min = 2.5xlO~ 5 G, and U B =~ 2.4 x 10~ 11 erg cm -3 . 

Considering the radial dependence of the optical photon field, we obtain a maximum 
predicted Inverse Compton emission due to scattering of the starlight ranging from 3.7 x 
10 36 erg s _1 at 1 pc galactocentric radius, to 7.9 x 10 35 erg s^ 1 at 500 pc (the maximum jet 
extension). We conclude that Inverse Compton radiation would be a small contribution to 
the X-ray emission, suggesting that synchrotron may be the dominant emission mechanism 
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if the S-shape feature is indeed a jet. 

3.2. Is hot gas responsible for the hard emission? 

We examine here the possibility that the origin of the hard S-shaped emission is thermal. 
We consider two scenarios suggested by analogies with other galaxies: thermal emission from 
the jet-ISM interaction, as in NGC 1052 (Kadler et al 2004), and shocks in the ISM, resulting 
from an outburst of nuclear activity, as suggested in the case of NGC 4636 (Jones et al 2002). 

An elongated feature in the central galaxy region, of temperature kT ~0.5 keV, was 
found with Chandra in NGC 1052, the prototype elliptical galaxy LINER (also at a distance 
of 22.6 Mpc, Knapp et al. 1978). This feature has 10 times the flux of the S-shaped emission 
in NGC 821 (~3.5xl0~ 13 erg cm^s^ 1 ), and about double the linear size (Kadler et al. 2004). 
In NGC 1052 the extended emission has radio and optical counterparts, and may be more 
fan shaped than linear, though the photon statistics are limited in this short (2.3 ksec) ACIS 
observation. Kadler et al. suggest shock heating of gas, converting some of the kinetic energy 
of the observed radio jet into X-ray emission. In the case of NGC 821, if a jet is present, 
it could similarly deposit energy in the surrounding ISM, causing it to produce the hard 
thermal emission. However, if its kinetic power is of the order of that estimated in Sect. 3.1 
(~ 10 37 erg s _1 ), this mechanism cannot explain the luminosity of the hard thermal emission 
measured with ACIS (L x ~ 1.9 x 10 39 erg s _1 , Table 2). Unless a jet is present with orders 
of magnitude larger kinetic power, other sources of energy deposition are needed. 

An alternative analogy is offered by the hot 'arms' of NGC 4636 (Jones et al 2002), 
a giant elliptical in Virgo with no reported nuclear activity or jets. These arms, having a 
larger spatial scale (8 kpc) and a lower temperature (~ 0.5 — 0.7 keV) than the S-shaped 
feature of NGC 821, are two symmetric features crossing the galaxy center, discovered in the 
Chandra ACIS data of NGC 4636. The NGC 4636 arms are accompanied by a temperature 
increase with respect to the surrounding hot ISM, which led to the suggestion by Jones et 
al of shock heating of the ISM caused by a nuclear outburst. We will discuss this process at 
the end of Sect. 3.3 below, in the unsteady accretion hypothesis. As discussed in § 2.3, the 
S-shaped feature of NGC 821 could similarly be hotter than its surroundings. Assuming a 
temperature of kT=3 keV for this feature, close to the lower limit suggested by our spectral 
fit (Table 2), we obtain a density n = 8. 59^};^ cur 3 (errors at 90%). Using the best-fit value 
(and a temperature of 3 keV), we obtain a thermal pressure exceeding by a factor of ~ 14 
that quoted in §3.1 for the surrounding hot ISM, suggesting a non-equilibrium situation, 
if only thermal pressures are involved. However, given the uncertainties in both T and n, 
the thermal pressures of both S-shaped feature and surrounding ISM could be similar. It is 
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clear that deeper Chandra observations are needed to better constrain the energetics of this 
feature. 



3.3. Is accretion present? 

Whether it is a jet or hot shocked gas, the S-shaped feature of NGC 821 requires a 
considerable energy input, an obvious source of which is accretion onto the nuclear SMBH. 
Taking at face value the indication of the two-component fit of the circum-nuclear diffuse 
emission, which is consistent with the presence of a kT ~ 0.5 keV thermal component, we 
can estimate the nuclear accretion rate, in the steady spherical accretion scenario of Bondi 
(1952). This estimate is rough, because the gravitational capture radius (which depends 
on the gas temperature and on the mass of the SMBH; see the textbook by Frank, King 
& Raine 2002) is r acc ~ 3 — 23 pc in our case, smaller than the physical resolution of the 
image (55 pc, Table 1). In this estimate of r acc we have allowed for the uncertainties in 
both the SMBH mass and kT. Based on the circum-nuclear gas temperature and density, 
that we estimate from the emission measure of the diffuse thermal component to be n = 
4.1+} 5 9 x 10~ 3 cm~ 3 , the Bondi mass accretion rate (also following Frank, King & Raine 2002) 
is M aC c = 1.1 x 10~ 7 — 2.0 x 10~ 5 M yr -1 , including again all the uncertainties in the SMBH 
mass, kT and n. The corresponding luminosity is L acc = 6.2 x 10 38 — 1.1 x 10 41 ergs s _1 , 
with the customary assumption of 10% accretion efficiency The luminosity of the S-shaped 
feature is within this range, therefore in principle it could be explained by Bondi accretion 
of the hot ISM. This conclusion, however, is most certainly not correct, because comparison 
with independent optical data, discussed below, shows that our uncertain X-ray data may 
have led to underestimating the circum-nuclear gas density and therefore M acc . 

Since the hot ISM is the thermalized integrated result of the stellar mass losses, as a 
minimum one expects the total stellar mass loss rate within r acc to be accreted (there 
may also be gas inflowing from outside r acc ). M* can easily be obtained from the luminosity 
density profile recovered from HST data for the central galaxy region (Gebhardt et al. 2003). 
Using a conversion factor from luminosity to mass loss rate for an old stellar population at the 
present epoch (e.g., Ciotti et al. 1991), this leads to M* = 9.8 x 10~ 6 and 2.6 x 10~ 4 Moyr" 1 , 
for the two extreme values of r acc . These M* values are respectively a factor of ~90 and 
~13 larger than the M acc values derived above for the same r acc . Given that this estimate of 
M* is quite robust, we must conclude that either the derivation of M acc above is inaccurate, 
or the gas is not steadily inflowing within r acc . The former possibility cannot be excluded 
with the present data, since in our calculation of M acc we used a density n value that is an 
average measured over a region extending much farther out than r acc ; n(r acc ) is likely to be 
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significantly higher than this average (e.g., a factor of ~ 30 times higher, for a n oc r °- 9 
profile) . 

Assuming that it is just the stellar mass loss rate within r acc that is steadily accreted, 
accretion luminosities > 20 times larger than the observed L x of the hard emission are recov- 
ered (from L acc = 0.1M*c 2 ). Then, the "solutions" proposed for the other X-ray faint nearby 
nuclei must be revisited for NGC 821. These are (a) accretion occurs but with low radiative 
efficiency (e.g., Narayan 2002); (b) accretion sustains a jet (this can be coupled again to a 
low radiative efficiency), whose total power is of the order of L acc , as in the modeling for 
IC 1459 (Fabbiano et al. 2003) and M87 (Di Matteo et al. 2003); (c) accretion is unsteady 
and therefore the hot ISM in the nuclear region needs not be inflowing (Siemiginowska et 
al 1996; Janiuk et al 2004). In this case the feedback from the central SMBH can be either 
radiative (Ciotti & Ostriker 2001) or mechanical (Binney & Tabor 1995, Omma et al. 2004) 
and make accretion undergo activity cycles: while active, the central engine heats the sur- 
rounding ISM, so that radiative cooling - and accretion - are offset; then the central engine 
turns off, until the ISM starts cooling again and accretion resumes. NGC 821 may be in a 
stage of such a cycle when a nuclear outburst has recently occurred. Note that the accretion 
luminosity that is radiatively absorbed by the ISM during an outburst (Ciotti & Ostriker 
2001) largely exceeds the hard thermal emission observed at the center of NGC 821. 

The presence of the central S-shaped feature that is hotter than the surrounding gas is 
uncontroversial evidence that central heating is at work, and therefore some type of feedback 
from the SMBH is occurring. The unsteady scenario seems promising to fit adequately the 
case of NGC821, also because this galaxy is clearly hot gas poor, as if recently swept by an 
outburst-driven wind. From our ACIS-S data we estimate un upper limit of ~ 4 x 10 6 M Q on 
the amount of hot ISM, many orders of magnitude smaller than for hot gas rich ellipticals 
(see Fabbiano 1989). 

4. Conclusions 

We have shown that NGC 821, which is used as template of a quiescent and normal 'old 
elliptical' (see, e.g. Ho et al 2003) shows clear signs of energy deposition by the nucleus: an 
elongated, bent, kiloparsec-size S-shaped feature, centered (within the ~ 1.5" errors) on the 
nucleus. This feature has a (0.3-10 keV) L x ~ 1-9 — 2.6 x 10 39 ergs s -1 , and a hard spectrum 
consistent with a power-law with T ~ 1.8, or thermal emission with kT > 2 keV. It may be 
embedded in a faint hot ISM with kT ~ 0.5 keV. 

This feature could represent the faintest yet reported extragalactic jet. We can exclude 
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Inverse Compton (either of radio photons or of the central stellar photon field) as an explana- 
tion of the X-ray emission. If this feature is indeed a non-thermal jet, this leaves synchrotron 
as a possible explaination. The a ra( n -x < 0.7 would be consistent with this hypothesis. In 
the jet scenario, we constrain the activity timescale to be ~ 4 x 10 4 years. We also estimate 
that in this jet most of the power may be in X-ray radiation, rather than kinetic energy. 

Alternatively, the emission could be thermal. It is unlikely that the detected luminosity 
could result from interaction of a jet with the surrounding ISM in this galaxy (as may occur 
e.g., in NGC 1052, Kadler et al 2004) because the radio upper limits imply a total jet kinetic 
power orders of magnitude below the detected X-ray luminosity. A more likely scenario 
is that of S-shaped shocks of the hot ISM, driven by an outburst of nuclear activity, as 
suggested in the case of NGC 4636 (Jones et al 2002). The characteristics of this feature 
(luminosity and hard spectrum) suggest that it is driven by nuclear activity. Accretion 
from stellar mass loss rate in the circum-nuclear region would produce a luminosity well in 
excess of the detected one, suggesting either low radiative efficiency (e.g., Narayan 2002) or 
unsteady accretion (e.g., Binney & Tabor 1995; Siemiginowska et al 1996; Ciotti & Ostriker 
2001). The second scenario is supported by the very small amount of hot ISM that may be 
present in NGC 821. It remains unexplored by the current modeling of outbursts through 
hydrodynamical simulations whether and how an S-shaped structure can be created. 

It is only in the X-ray band, thanks to the Chandra spatial resolution, that signs of 
nuclear activity are found in NGC 821. Either jet or shock explanation of the S-shaped 
nuclear feature detected with Chandra suggest a fairly recent nuclear outburst, now spent. 
However, the present data are not deep enough to either allow a detailed study of the 
morphology and spectral parameters of this feature, or a model-independent detection of the 
circum-nuclear hot ISM. A significantly deeper Chandra exposure will be needed to answer 
the many questions raised by the present data. Deeper VLA data are also needed to set a 
tighter constraint to the jet scenario. 

We thank the Chandra X-ray Center DS and SDS teams for their efforts in reducing the 
data and developing the software used for the data reduction (SDP) and analysis (CIAO). We 
have used the NASA funded services NED and ADS, and browsed the Hubble archive. This 
work was supported by NASA contract NAS 8-39073 (CXC) and NASA grant G03-4133X. 
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Fig. 1. — True color adaptively smoothed image of NGC 821. The ellipse marks the 25th 
mag isophote. 
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Fig. 2. — a) True color image of the central region of NGC 821, unsmoothed. The cross 
and surrounding circle represent the 2MASS nuclear position and uncertainty, from NED. 
b) A larger field image showing both the wavdetect source regions (yellow) and the spectral 
extraction regions for the 'jet' and the diffuse emission (light blue). 
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HR1 ([0.5-1] keV vs. [1-2] keV) 

Fig. 3. — X-ray colour-colour diagram of 'jet' and diffuse emission regions. Typical hardness 
ratios for a power-law model (r = 1 — 2.5) are plotted in blue. Typical hardness ratios for 
a thermal model (kT = 0.25 — 20 keV) are plotted in red. Galactic line-of-sight absorption 
(Nh = 6.4 x 10 20 cm -2 , Stark et al 1992) is assumed in both models. 
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